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Abstract 



Motivated by the problem of coupled Hubbard chains, we investigate a generalisation 
of the Schulz-Shastry model containing two species of one-dimensional fermions interacting 
via a gauge field that depends on the positions of all the particles of the other species. 
The exact many body ground state of the model can be easily obtained through a unitary 
transformation of the model. The correlation functions are Luttinger-like - i.e., they decay 
through power laws with non-integer exponents. Through the interaction dependent cor- 
relation functions of the two-particle operators, we identify the relevant perturbations and 
hence, possible instabilities. 
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Exactly solvable models |1|, 0, 0] have always attracted a lot of interest in theoretical 
physics, because they serve as paradigms for more complicated systems. The fact that 
these models are usually in one dimension no longer make them unrealistic, however, since 
current technological advances have seen the advent of many semi- artificial one dimensional 
systems, such as quantum wires, quantum Hall bars, one- dimensional organic metals and one 
dimensional spin chains 0. In fact, phenomena such as one dimensional Luttinger liquids and 
the Haldane gap in spin chain models have actually been experimentally seen]^. Besides their 
role in these systems, exactly solvable models have played a very important role as a reliable 
test for various approximation methods and for developing qualitative understandingP]. 

However, for two or more dimensions, there have been very few exact results. For in- 
stance, the large U Hubbard model has been studied using several approximation schemes|0], 
none of which have led to completely reliable results. In recent years, there have been at- 
tempts to understand two dimensions through the coupling of one-dimensional chains. Both 
coupled spin chains]^ and coupled Hubbard models with interchain hopping and inter- 
chain interactions have been studied using a variety different schemes such as weak coupling 
renomalisation group techniques and bosonisationfl^, exact numerical diagonalisations [|Tl| , 



etc. Unfortunately, in the absence of any exact results, the interpretation of the results of 
these inter-chain coupling studies have remained difficult fT^ . 



With the motivation of approaching two dimensional phenomena through the coupling 
of one- dimensional chains, in this paper, we study a generalisation of a class of models [|13| 
that can be diagonalised by a pseudo-unitary transformation and still exhibit non-trivial 
Luttinger liquid behaviour. Our model has two species of particles with pseudospin index 
(T = ±, at the positions x^i and with momenta Pai, with a Hamiltonian given by 



N 

'-ai ) 



Here, Hg-j = p^-j + crAo-(xo-j) is the 'covariant momentum' introduced in Ref.|]T3[ and is a 
'chain index' or 'band index'. (The nomenclature of chain or band index will be explained 
later.) We have chosen to have only even powers of the covariant momentum in the Hamil- 
tonian, although positivity of energy only requires that the largest power of the covariant 
momentum be even. This maintains the symmetry x — > —x or parity, which simplifies the 
presentation of the calculations, although the result goes through even when we include odd 
powers of the momenta. As explained in Ref. ||13||, particles interact via a gauge potential, 
given for the particle at the position x by A^ji^x) = Y^j V{x — X-^j) - i.e., the potential for 
the particles with positive pseudospin is due to the presence of the particles with negative 
pseudospin and vice-versa. The potential is chosen to be an even function, vanishes at infin- 
ity and explicitly breaks time-reversal invar iance, although it is invariant under a combined 
operation of time reversal and reversal of pseudo-spin index. We have generalised the model 
in Ref. [0 by including a chain index (or equivalently a band index) and allowing higher 
powers of the covariant momentum in the Hamiltonian. Our model reduces to the Schulz- 
Shastry model for J = 1 and ai = 1. Clearly, the aj are not dimensionless, and in fact, 
explicitly contain a scale A (except for oi, which is dimensionless). 
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As noted by the authors in Ref.[jT3[, the same pseudo-unitary transformation that they 
use to diagonahse their Hamiltonian, 



diagonahses any power of Uai, as long as we choose the function S (a function of the 2n posi- 
tions of the particles) to eliminate the interaction in Eq.(p. Thus, we obtain the transformed 
Hamiltonian given by 

N 

jj ^ gi5({x+.},{x_j)^g-i5({x+j,{x_j) =Y^J2 aiip^if = H^, (3) 

where the interaction pieces have been removed by the transformation. However, the eigen- 
values and eigenfunctions are not the same as that for a genuinely non-interacting Hamilto- 
nian because the boundary conditions on the wave-functions are now different. 

For the single particle Hamiltonian if^j in Eq.@, the eigenvalue equation is a 2A^*^ order 
differential equation, and depending on the energy chosen, will have at most 2N different 
solutions. The general solution is given by 

N 

= Y,cie"'"' + h.c., (4) 
I 

where the fc/'s are known in terms of E and the A^ — 1 constants aj (we always choose a^r = 1 
without loss of generality since it only sets the overall scale). However, not all fc/'s need be 
integer multiples of 27r/L, where L is the size of the system. For those that are not, the 
corresponding c/ vanish so as to make the wave function periodic in L. Since the a/'s are 
fixed, we may choose only one of the fc/'s to be independent, say k, which in turn fixes the 
dispersion to be 

N 

E{k) = Y.ai{kY . (5) 
I 

The Fermi points are the roots of the equation E{k) = Ep- We choose an energy Ep where 
the 2A^ roots {—bj, 6/} are all real and distinct (with bi < b2 < ■ ■ ■ < &Ar). The Fermi points 
{—Ei,Ej} are given by Ej = 27T[ri]/L where bj = 27irj/L and [rj] stands for the largest 
integer below rj. 

We now see the justification for calling I the chain or band index. The usual identification 
of an iV-chain model with an iV-band model is made by diagonalising the kinetic energy of 
the A^-chain model and using the A^ different momenta along the direction perpendicular to 



the chain to label the A^ bands or A^ dispersion relationsfT^. The filling of these bands upto 
the Fermi level defines the set of 2A^ Fermi points ±kp. Our model is slightly different from 
the A^-chain model in that it has only a single dispersion(Eq.(^). However, the dispersion 
is not quadratic and has A^ wells (unlike the usual quadratic dispersion which has one well 
per band) and 2A^ Fermi points. Hence, the physics it describes is similar to that of the A^ 
band model. 
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We are interested in the solution of the original Hamiltonian in Eq.(|l]) and not the 
transformed Hamiltonian in Eq.(^). Although the single particle energies of the two Hamil- 
tonians are the same, their wave-functions are related by the pseudo-unitary transformation 
if) = e~'^^ip, where S was chosen to cancel the interaction and is of the form 

S{{x+i}, {x-i}) = E{x+i ~ x^j) where E{x) = / dx'V{x'). (6) 

We can compute the difference between S{x_i = L) and S{x_i = 0) for any particular 
negative pseudospin coordinate 

S{x_, = L)-S{x_, = 0) = Y.[E{x+j - L) - E{x+,)] (7) 

j 

= nl V{x)dx = nl5. (8) 

in terms of a phase shift S and which is the total number of positive pseudospin particles. 
One gets a similar result if we choose the reference particle to be a positive pseudospin 
particle, with the only difference that gets replaced by and 6 by —6. Hence the 
quantisation condition on the wave-numbers of the particles becomes 

Lk±i =1= n^S = 27cn±i (9) 

where the n±i are integer quantum numbers analogous to those used in the non-interacting 
case. Since, in general, ri^S ^ integral multiple of 27r, the free Hamiltonian and the inter- 
acting Hamiltonian are in different Hilbert spaces. 

So far, all the arguments used by Schulz and Shastry have gone through for our model 
as well. The differences begin when we try to construct the many body ground state and 
the spectrum of low energy excitations. For ease of presentation, we will now specialise to 
the two band case, explicitly perform the calculations leading to the low energy effective 
Hamiltonian, and then generalise to the case of bands. 

For the two-band case, the single particle dispersion is given by 

E = e^'He-^' = ^ a,pt, + a,pl, (10) 

where 04 = has length dimension two and 02 = —1 is dimensionless. As mentioned before, 
we restrict the Fermi level to lie within the double well - i.e., we have four distinct Fermi 
points (— F2, — Fi) on the left and (^1,^2) on the right. For each energy, the degeneracy 
is either four or two , depending on whether or not both ki and ^2 satisfies the boundary 
condition given in Eq.(^). However, all that really matters is that energy levels in both the 
wells below the Fermi level are filled. Let us assume that in band 1, there are rii states below 
the Fermi level and in band 2, there are n2 states below the Fermi level for both pseudospins. 
The ground state energy is then given by 

e,s=2y: e[m'-^2^'] (11) 

o-=± i=— m 
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where hi = A^(^)^ and h2 = (x)^- have assumed that n^5/27r is an integer and the 
factor of two takes care of the contributions from both the right and left moving sectors. 
(Note that the states on the right branches of both the wells are right movers, whereas the 
states on the left branches of both wells are left movers. See figure below.) 




k=0 



• Fig 1. Dispersion for the two band model. For the Fermi energy Ep, the two bands 
are denoted by dotted indices (band 1) and full lines (band 2). {—F2,Fi) are clearly 
left-mover Fermi points and (—^1,^2) are right mover Fermi points. The inset shows 
the dispersion of the three band model, with left-mover branches denoted by full fines 
and right-mover branches denoted by dotted lines. 

The second order fluctuation in the energy due to the addition of nji±i and nL±i particles 
for the right and left movers can also be computed. It is given by 

{n2+nR±2) 

E^'^ = [E E Mt + c^r-h^it + c^n+R^L], (12) 

± -(ni-nK±i) 
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= EE [hifiM + hMnj)] [riR^i-c^f (13) 
= EE^?K) [J±i±n^S/2nr + in^jr] (14) 

where = 5/271 in the first row and the subscript 2 is to indicate that we keep only terms 
upto quadratic order in the fluctuations. In the second row, /i(n/) = n/(n/ + l)(2r;,/ + 1) 
and /2 = ri/ + 1/2, and in the third row, we have defined the current J±j = nji±j — nL±j, the 
charge n±j = nji±i + ni^i and the 'density' g{ni) = hifi{ni) + h2f2{ni). The total charge 
is clearly = rij-i + n±2- Note however, that unlike the Schulz-Shastry model, here fourth 
order fluctuations do exist, which we neglect because we are only interested in low energy 
fluctuations. 

We bosonise as in the single band case by introducing boson fields (j)±i with their conjugate 
momenta Il±i. These are related to the currents and charge densities as 

n±i = -^d^(f)±i (15) 

'71 



J±i = -^T^±I. (16) 

'71 



(We use the notation of Ref. [[L4||.) To rewrite the effective hamiltonian for the low energy 
fluctuations in terms of the boson fields, we have to identify the function g{nj)L/27i = pj as 
an effective density after which we obtain 

H = jzT.I dxpi{[-I[i^ ± -{dAiT + ^x02t]' + {dA±i?}. (17) 

But interestingly, although pi contains information about the scale, the low energy effective 
Hamiltonian is scale invariant - there are no mass terms (or cosine terms leading to mass 
terms) for the boson fields. A similar redefinition of variables as in the one band case, 

<i>±I = <t>±I , = T -{d^<P^i + da,<p^2) (18) 

71 

leads to a non-interacting form of the Hamiltonian given by 

H = j2Y.fd^pi [(n±/)' + {dj±if ]■ (19) 
I ± 

Thus the correlators of the tilde fields are just free field correlators. In terms of the non- 
tilde bosonic variables or equivalently in terms of the fermion fields, the Hamiltonian is not 
non-interacting. However, since they are explicitly known in terms of the free fields, their 
correlators can also be explicitly calculated. 

In fact, at this stage, the generalisation to N chains is obvious. The single particle 
dispersion of the iV-band model has wells and 2A^ Fermi points. The Hamiltonian for 
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quadratic fluctuations about the Fermi points is precisely the same as that in Eq.(|T7|) with 
the replacement 

N 

(<9x0Tl + dx(l)T2) J2 9x<PtJ (20) 
J 

As before, the redefinition of 0/ and 11/ in terms of the tilde fields leads to the non-interacting 
form of the Hamiltonian in Eq.(|T9|) with the sum going over all bands. 

We now compute correlation functions using the representation for the fermion operators 
in terms of the non-interacting boson fields given by 
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N 



i}R±i = exp(0ij±/ ^ — ^0^), (21) 
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N 



^L±i = exp(0i±,±— ^0^). (22) 



(We follow the notation in Ref. |Tj] and define 4>r±i = l/2{(f)±j — Il±i{x)dx) and (f)L±i = 
l/2(0±/ + n-|-/(x)(ix). ) The one-particle correlation function is given by 

Gz±i{x) =< fe/(x)4±,(0) >~ x-^ (23) 

with 7] = 1 + N{6'^/2tt'^) for both right and left movers {Z = R/L), for both pseudospins 
and for all /. As in the one-band case, the fermion has an anomalous dimension given by 
T) 7^ integer. This is the indication that the system is a Luttinger liquid and not a Fermi 
liquid. The interesting point to note here is the dependence of the anomalous dimension 
on the number of chains. The model is not just a collection of one-band Luttinger liquids - 
there exists a genuine dependence on the number of bands. 

We can also compute exponents of the two-particle operators(TPO) in order to identify 
the relevant perturbations and hence potential singularities. In the one-band case, the only 
non-trivial two-particle correlations involved excitations at both the right and left Fermi 
points, because these were the only two Fermi points. Here, however, we can have non- 
trivial two-particle correlations involving excitations at two right Fermi points and two left 
Fermi points as well. These exponents for the two particle correlations are tabulated below. 



TWO PARTICLE CORRELATIONS 



TPO 


V 


TPO(J ^ J) 


V 




2 




2{1 + N^) 








2(1 + 1^) 


'ipR±l1pL^:J 


2(1 + Ar5) 


ipz±i'ipz^j 


2 


'i/jR±li^LTJ 


2(l±A) + iVfi 


'ipZ±l'ipZ^:J 


2(1 + f^) 



Fortunately, they are independent of the band index and only depend on whether they 
involve both right and left Fermi points or right (left) movers at both Fermi points. Interest- 
ingly, none of the RR or LL exponents lead to relevant perturbations. This is in agreement 
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with the weak couphng RG approach[0, where there is a non-zero contribution to the four- 
point vertex only when there is momentum transfer between left and right Fermi points. In 
the single chain case considered in Ref . , there was no possibility of momentum transfers 



between two Fermi points on the left or two on the right, since the model only had one on 
each side. However, even in the more general case of left-moving Fermi points and 
right moving Fermi points |jT6|, graphs involving loop momenta in two left-moving shells or 
two right moving shells are always zero, because the energies have the same sign and the 
contour integral for the energy vanishes. 

For positive 6, the only relevant perturbations are 

i^Ui^L-j // = 1 + (1 - + {N- l)5Vv^^ 

iPR-ii^L+j r] = l + {l- 5/Tif + {N- l)5Vvr'. (24) 

Clearly as the number of chains increases, the exponent increases until at some critical value 
of = A'^c oc 1/5, both the operators above cease to be relevant. Similarly, for negative 5, 
the relevant perturbations are 

^..V'L+J r/ = 1 + (1 + + (AT - 1)5^ 

iPr+i^l-j r/ = 1 + (1 + 6/71)' + {N- l)(5Vvr^ (25) 

which cease to be relevant beyond Nc- Hence, for N > Nc chains, there are no relevant 
perturbations at all. The system is always a Luttinger liquid. For A^ < Nc, the system, 
depending on which instability grows, (which perturbation will be added), which in turn, 
will be dictated by the realistic model that we wish to study, will be in different ground 
states. 

Let us compare our results with the results obtained by giving additional internal degrees 
of freedom to the o" = ± particles [1^. If we assume that they occur in m-flavours, then the 
Hamiltonian is just 

m 

^ = EE«Kn/..)'. (26) 

1=1 ai 

Surprisingly, an analogous calculation leads precisely to the same exponents as in Eqs.(^) 
and (|25|) with A^ replaced by m. However, in this case, there genuinely exist 2m degrees of 
freedom, and the various two particle correlators have physical meaning. I = J give two par- 
ticle correlators of the same particle, whereas there is no analogue of I ^ J correlators. The 
exponents are actually independent of the particle index because of the internal symmetry. 
For our Hamiltonian in Eq.(l), however, there are only two degrees of freedom correpond- 
ing to the a = ± particles. It is only after linearising around the different Fermi points 
and assuming that each of the linearised fermions can be bosonised independently that we 
have A^ independent R/L moving fermions or bosons, whose correlators can be computed 
independently. For the original fermions, the only relevant charges are and the relevant 
currents are J± = Y^i{nR±i — nL±i)- 

Many of the issues in coupled chain models, however, remain unaddressed in this rather 
simple model, which is perhaps better thought of as a single chain model with a more 
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complicated band structure. To really apply this model to A^-cliains, one would have to 
modify the model, so that there is some analog of the interplay between interchain hopping 
and intra-chain interactions. However, note that even as an A^-band model, it is not trivial 
that the correlation functions are identical to those of the m- flavour model. 

In conclusion, we have studied a general model with 2N Fermi points, (an A^-band or N 
chain model), which is exactly solvable and has non-trivial Luttinger liquid behaviour. We 
computed the exponents of the various two-particle operators and found the possible relevant 
perturbations. Interestingly, we found that the exponents have non-trivial dependence on 
the number of chains - they are not merely additive. Furthermore, beyond a certain number 
of chains iVc, which is inversely proportional to the strength of the interaction 5, all pertur- 
bations are irrelevant and the Luttinger liquid ground state is robust. Thus, this model, if it 
can be effectively generalised to higher dimensions would be a good starting point to study 
possible Luttinger liquid ground states in higher dimensions. 



Acknowledgments 

RKG would like to thank the Mehta Research Institute for hospitality during the course of 
this work. 



References 

H. A. Bethe, Z. Phys. 71,205 (1931); for a review and further references, see references 
i and i. 

J. M. Luttinger, J. Math. Phys. 4, 1154 (1963); F. D. M. Haldane, J. Phys. C 14,2585 
(1981); ibid, Phys. Rev. Lett.45, 1358 (1980). 

D. C. Mattis and E. H. Lieb, J. Math. Phys. 6, 304 (1965). 

For a recent introduction to one-dimensional systems, see Zachary N. C. Ha, 'Quantum 
Many-Body systems in One Dimension', World Scientific, Singapore, 1998. 

M. Azuma et al, Phys. Rev. Lett. 73, 3463 (1994); T. M. Rice and E. Dagotto, Science. 

V. J. Emery, in ' Correlated Electron Systems', ed. by V. J. Emery, World Scientific, 
Singapore, 1992. 

A. Houghton and J. B. Marston, Phys. Rev. B50, 1351 (1994); A. H. Castro-Neto and 

E. Fradkin, Phys. Rev. Lett. 72, 1393 (1994), J. R. Engelbrecht and M. Randeria, Phys. 
Rev. Lett. 65, 1032 (1990); M. Fabrizio et al, Phys. Rev. B44, 1033 (1991). 

H. J. Schulz, Phys. Rev. Lett. 77, 2790 (1996); ibid, Phys. Rev. B53, R2959, (1996); L. 
Bartosch and P. Kopietz, |cond-mat/9703205| . 



9 



[9] M. Fabrizio et al, Phys. Rev. B46, 3159 (1992); A. M. Finkelstein and A. R. Larkin, 
Phys. Rev. B47, 10461 (1993); A. Nersesyan, A. Luther and F. Kusmartsev, Phys. Lett. 
A 176, 363 (1993). 

[10] L. Balents and M. P. A. Fisher, Phys. Rev. B53, 121333 (1996); H. Lin, L. Balents and 
M. P. A. Fisher, |cond-mat/9703055| . 

[11] C. A. Hayward and D. Poilblanc, Phys. Rev. B53, 11721 (1996); H. Enders et al, Phys. 
Rev. B53, 5530 (1996); R. M. Noack et al, Europhys. Lett. 30, 163 (1995). 

[12] D. G. Clarke, S. P. Strong and P. W. Anderson, Phys. Rev. Lett. 74, 4499 (1995). 

[13] H. J. Schulz and B. S. Shastry, Phys. Rev. Lett. 80, 1924 (1998). 

[14] R.Shankar, 'Bosonisation: How to make it work for you in condensed matter' , in Low- 
Dimensional Quantum Field Theories for Condensed Matter Physics', World Scientific, 
1992. 

[15] R.Shankar, Rev. Mod. Phys. 66, 129 (1994). 
[16] P. K. Mohanty and S. Rao, work in progress. 



10 



